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ABSTRACT 

The morphology and kinematics of molecular clouds (MCs) are best explained as 
the consequence of super-sonic turbulence. Super-sonic turbulence fragments MCs into 
dense sheets, filaments and cores and large low density "voids" , via the action of highly 
radiative shocks. We refer to this process as turbulent fragmentation. 

In this work we derive the mass distribution of gravitationally unstable cores gen- 
erated by the process of turbulent fragmentation. The mass distribution above one 
solar mass depends primarily on the power spectrum of the turbulent flow and on the 
jump conditions for isothermal shocks in a magnetized gas. For a power spectrum index 
j3 = 1.74, consistent with Larson's velocity dispersion-size relation as well as with new 
numerical and analytic results on super-sonic turbulence, we obtain a power law mass 
distribution of dense cores with a slope equal to 3/(4 — /3) = 1.33, consistent with the 
slope of the stellar IMF. Below one solar mass, the mass distribution flattens and turns 
around at a fraction of a solar mass, as observed for the stellar IMF in a number of 
stellar clusters, because only the densest cores are gravitationally unstable. The mass 
distribution at low masses is determined by the probability distribution of the gas den- 
sity, which is known to be approximately Log-Normal for an isothermal turbulent gas. 
The intermittent nature of the turbulent density distribution is thus responsible for the 
existence of a significant number of small collapsing cores, even of sub-stellar mass. 

Since turbulent fragmentation is unavoidable in super-sonically turbulent molecular 
clouds, and given the success of the present model in predicting the observed shape of 
the stellar IMF, we conclude that turbulent fragmentation is essential to the origin of 
the stellar IMF. 

Subject headings: turbulence - ISM: kinematics and dynamics - stars: formation - 
stars: mass function 
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1. Introduction 

The process of star formation, particularly the ori- 
gin of the stellar initial mass function (IMF), is a 
fundamental problem in astrophysics. Photometric 
properties and chemical evolution of galaxies depend 
on their stellar content. The process of galaxy forma- 
tion cannot be described independently of the process 
of star formation, since galaxies are partly made of 
stars. Stars are also an important energy source for 
the interstellar medium of galaxies. 

Stars are formed in molecular clouds (MCs), which 
have been the focus of the research on star formation 
for more than two decades. Currently, there is no 
generally accepted theory of star formation, capable 
of predicting the star formation rate and the stellar 
IMF based on the physical properties of MCs. This is 
hardly surprising, since turbulent motions are ubiqui- 
tously observed in MCs, and the physics of turbulence 
is poorly understood, due to the great mathemati- 
cal complexity of the fluid equations. Magnetic field, 
self-gravity and high Mach numbers further increase 
the complexity. 

The steady growth of computer performance has 
now made large three— dimensional numerical simula- 
tions of super sonic magneto-hydrodynamic (MHD) 
turbulence feasible (Padoan & Nordlund 1997, 1999; 
Stone, Ostriker & Gammie 1998; Mac Low et al. 1998; 
Padoan, Zweibel & Nordlund 2000; Klessen, Heitsch 
& Mac Low 2000; Mac Low & Ossenkopf 2000; Os- 
triker, Stone & Gammie 2000; Heitsch, Mac Low & 
Klessen 2000; Padoan et al. 2001a, b). Comparisons of 
numerical experiments with observational data have 
shown that super-sonic turbulence can explain the 
morphology and kinematics of MCs, and the forma- 
tion of dense cores, provided that the motions are 
also super- Alfvenic (Padoan, Jones & Nordlund 1997; 
Padoan et al. 1998; Padoan & Nordlund 1997, 1999; 
Padoan et al. 1999; Padoan, Rosolowsky & Goodman 
2001; Padoan et al. 2001a, b). We refer to this process 
of formation of dense cores in MCs by super sonic tur- 
bulence as turbulent fragmentation. Since protostars 
evolve from the collapse of gravitationally unstable 
cores in MCs, even the stellar IMF could then be the 
result of turbulent fragmentation, with the power law 
shape of the IMF ultimately being the consequence of 
the self-similar nature of turbulence. 

In this paper we do not use these increasingly 
sophisticated numerical simulations directly, but in- 
stead develop an analytic model. The assumptions of 



the model are inspired by the qualitative^ properties of 
the numerical model, but the current work does not 
depend on any particular set of numerical models or 
results. 

Previous analytic models by Larson (1992), Hen- 
riksen (1986, 1991) and Elmegreen (1993, 1997, 1999, 
2000b) have derived the stellar IMF on the basis of the 
self-similar structure of MCs. Larson, assuming one 
dimensional accretion, predicted a rather steep IMF 
slope, equal to the MC fractal dimension, while Hen- 
riksen found the IMF slope to depend on both the MC 
fractal dimension and the relation between density 
and linear size for structures inside MCs. Elmegreen 
pointed out that the IMF that results from random 
sampling of a self-similar cloud has an exponent x = 1 
(Salpeter x = 1.35), independent of the cloud fractal 
dimension. In Elmegreen's model the IMF is steeper 
than X = 1 because the random sampling rate is as- 
sumed to be proportional to the square root of den- 
sity, and because of "mass competition" . 

These models are based on assumptions about the 
cloud geometry justified by the apparent fractal struc- 
ture of MCs (Beech 1987; Bazell & Desert 1988, Scalo 
1990; Dickman, Horvath & Margulis 1990; Falgar- 
one, Phillips & Walker 1991; Zimmermann, Stutzki & 
Winnewisser 1992; Henriksen 1991; Hetem & Lepine 
1993; Vogelaar & Wakker 1994; Elmegreen & Falgar- 
one 1996), but the processes responsible for generat- 
ing the assiimed geometry are not discussed in detail. 

In a previous attempt to relate the stellar IMF 
to the physical properties of super-sonic turbulence 
(Padoan, Nordlund & Jones 1997), we obtained the 
distribution of the local Jeans' mass in super-sonic 
isothermal turbulence, from the probability distribu- 
tion of the gas density. That work also provides a 
prediction for the stellar IMF, by identifying each lo- 
cal Jeans' mass with a protostar. Although the pre- 
diction of the lower mass cutoff and the low mass 
portion of the IMF might be roughly correct, this 
model under-estimates the number of massive stars 
relative to low mass stars (as any other Log-Normal 
IMF). The main reason, as pointed out by Scalo et 
al. (1998), is the unphysical assumption that the most 
massive stars originate from gas at relatively low den- 
sity including only a small fraction of the total mass. 

In the present work we consider specific proper- 
ties of MC turbulence and their consequences for the 
formation of protostellar cores. In particular, we 
i) assume approximate self-similarity of the super- 
sonic and super-Alfvenic velocity field, expressed by 
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a power law shape of the power speetnim of turbu- 
lence, and ii) assume that the jump conditions for 
isothermal MHD shocks determine the scaling of the 
typical size of protostellar cores. Based on these 
fundamental assumptions, we are able to derive a 
power law mass distribution of gravitationally unsta- 
ble cores, with a slope consistent with the Salpeter 
stellar IMF (Salpeter 1955). In order to derive the 
mass distribution of gravitationally unstable cores be- 
low one solar mass we also make use of the Probabil- 
ity Density Function (PDF) of the mass density in 
super-sonic turbulence, following Padoan, Nordlund 
& Jones (1997). 

The formation of dense cores and the relation be- 
tween the size of cores and the thickness of the post- 
shock gas are discussed in the next section. In § 3 
the power spectrum of super-sonic turbulence is dis- 
cussed. Results on the PDF of the gas density in 
super-sonic turbulence are summarized in § 4. The 
mass distribution of dense cores formed by the pro- 
cess of turbulent fragmentation is derived in § 5, and 
the mass distribution of gravitationally unstable cores 
is computed in § 6. The relevance of these results for 
the origin of the stellar IMF is discussed in § 7 and 8. 
Conclusions are summarized in § 9. 

2. The Origin of Dense Cores in Super— Sonic 
Turbulence 

The qualitative properties of the local density max- 
ima that form in super-sonic turbulence follow from 
first principles and numerical experiments have con- 
firmed such properties. 

Dense cores are formed in roughly isothermal super- 
sonic turbulent flows as the densest parts of sheets or 
filaments of shocked gas. Their typical size is there- 
fore comparable to the sheet thickness, A. Assum- 
ing that the magnetic pressure in the postshock gas 
exceeds the thermal pressure (this is the case for a 
large range of values of the preshock magnetic field 
strength, even when the flow is super- Alfvenic), the 
isothermal shock jump conditions are; 
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where poi and pi, Bi are the values of the gas den- 
sity and magnetic field strength before and after the 
shocks, respectively. L is the linear extension of the 
gas before the shock (measured in the direction per- 
pendicular to the shock surface) and A is the thickness 
of the postshock gas. A^q is the Alfvenic Mach num- 
ber of the shock, that is the ratio of the flow velocity 
and the Alfven velocity measured in the preshock gas: 
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(4) 



Here the relevant magnetic field components are the 
ones parallel to the shock surface, because the per- 
pendicular component does not provide pressure sup- 
port against the compression and it is not amplified 
by the compression. Since only the parallel compo- 
nents of B are amplified, the field in the postshock 
sheets is nearly parallel to the sheet surface, and elon- 
gated in the direction of dense filaments (real ones, or 
two-dimensional sections of sheets). Since dense cores 
are formed predominantly in corrugations of sheets, 
the magnetic field within dense cores can occasion- 
ally show a strong curvature, depending on the ori- 
entation relative to the line of sight. This may con- 
siderably complicate the interpretation of dust polar- 
ization measurements (Ward-Thompson et al. 2000) 
and should be taken into account. 

3. The Power Spectrum of Super— Sonic Tur- 
bulence 

The power spectrum of turbulence in the inertial 
range (below the energy injection scale and above the 
dissipation scale) may be assumed to be a power law. 
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where k is the wave number, and the spectral in- 
dex is /3 « 5/3 for incompressible turbulence (Kol- 
mogorov power spectrum) (Kolmogorov 1941), and 
/3 « 2 for pressureless turbulence (Burgers power 
spectrum) (Burgers 1974; Gotoh & Kraichnan 1993). 
The /3 = 2 power spectrum has often been assumed 
to be the correct power spectrum of super-sonic tur- 
bulence, at least for the compressional component of 
the velocity field. 

The purpose of the present work is primarily that 
of establishing analytically a relation between the 
power spectrum of turbulence and the stellar IMF 
and we only need to assume that the power spectrum 
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is a power law in the inortial range. Our analytic 
model does not assume a specific value of the spectral 
index, as obtained for example in numerical simula- 
tions. However, the predictions of our analytic model 
will be tested using the power spectrum of super-sonic 
turbulence computed from numerical simulations or 
estimated from observational data. For that reason 
we briefly summarize in this section some numerical 
and observational results. 

3.1. The Power Spectrum in Numerical Sim- 
ulations 

The most detailed study of the power spectrum 

of numerical compressible turbulence has been pre- 
sented by Porter, Pouquet & Woodward (1992, 1994) 
and Porter, Woodward & Pouquet (1998). From the 
1992 paper to the 1998 paper the largest numerical 
resolution increased from 256"^ to 1024"^. These works 
are limited to decaying turbulence, with Mach num- 
bers close to unity initially, and below unity at later 
times. The runs are therefore sub-sonic, except for 
an initial period of time. A magnetic field is not in- 
cluded. 

The velocity field is usually decomposed into its 
solenoidal v* and compressional v'^ components as 
V = V* + v= with V • = and V x v'' = 0. 
The velocity Fourier spectrum, E{k), is also sep- 
arated into its solenoidal and compressional parts: 
E{k)=E-'{k) + E%k). 

Porter, Pouquet & Woodward (1992) found that 
the compressional modes have a power spectrum 
E''{k) oc fc~^, and the solenoidal modes E^{k) oc k^^. 
In the later works, after the larger numerical resolu- 
tion runs were performed, the same authors concluded 
instead that both compressional and solenoidal modes 
develop a Kolmogorov power spectrum, E'^{k) oc 
E^{k) oc with E^/E^ « 0.15 (Porter, Pouquet 

& Woodward 1994; Porter, Woodward & Pouquet 
1998). This conclusion is not strongly supported by 
the plots of power spectra presented by the authors, 
since the largest resolution runs (512^ and 1024'^) are 
consistent with a Kolmogorov power spectrum for the 
solenoidal modes only over a very limited range of 
wave-numbers, approximately 4 < fc < 10. At larger 
wave-numbers, the power spectrum is flatter, approx- 
imately E'^{k) oc k~^. An interpretation of the shal- 
lower power spectrum of solenoidal modes at large 
wave numbers is provided, in terms of a "near dissi- 
pation range" . The discrepancy between the steeper 
(Burgers) power spectrum of compressional modes in 



the early 256'^ runs (Porter, Pouquet & Woodward 
1992) and the Kolmogorov power spectrum in the lat- 
est 256^ runs (Porter, Pouquet & Woodward 1994; 
Porter, Woodward & Pouquet 1998) is not discussed. 

We have recently started to perform a large num- 
ber of numerical experiments with the purpose of 
computing the spectral index of the inertial range of 
driven super-sonic MHD turbulence as a function of 
the sonic and Alfvenic rms Mach numbers of the flow. 
The experiments use an isothermal equation of state, 
uniform initial density and magnetic flelds, random 
initial velocity and random large scale forcing (both 
solenoidal), as described in previous works (Padoan et 
al. 1998; Padoan & Nordlund 1999; Padoan, Zweibel 
& Nordlund 2000). The results of these new experi- 
ments will be reported and discussed elsewhere. Re- 
sults on power spectra and structure functions are 
reported in Boldyrev, Nordlund & Padoan (2002), 
where the scaling relations are also predicted on the 
basis of a new analytic model of super-sonic turbu- 
lence (Boldyrev 2002). A power spectrum intermedi- 
ate between the Burgers and the Kolmogorov power 
spectra is found, E{k) oc A;~^'^^. 

3.2. The Power Spectrum from Observational 
Data 

The observed velocity dispersion-size Larson rela- 
tion (Larson 1979, 1981), Aw oc L", should reflect 
the power spectrum of turbulence in the interstellar 
medium (a = (/? — 1)/2), although it is often obtained 
from a combination of different clouds and cores inside 
the same cloud, or even different molecular transitions 
(see Goodman et al. 1998 for a discussion of the line 
width-size relation). Larson flnds a = 0.37 in the 
range of scales 1 < i < 1000 pc (Larson 1979); and 
a = 0.38 in the range of scales 0.1 < i < 100 pc 
(Larson 1981); Leung, Kutner & Mead (1982) ob- 
tain a = 0.48 for 0.2 < L < 4 pc; Myers (1983) 
gets a = 0.5 for 0.04 < i < 10 pc; Sanders, Scov- 
ille & Solomon (1985) find an unusually large value 
a = 0.62 for 20 < L < 100 pc, which they use to 
rule out any relation between a turbulent power spec- 
trum and the line width-size relation; Dame et al. 
(1986) obtain a = 0.5 for 10 < L < 150 pc; finally 
Falgarone, Puget & Perault (1992) use a compilation 
of data from the literature together with their own 
new data, in order to sample a very large range of 
scales, 0.01 < L < 100 pc, and include also a signifi- 
cant number of unbound objects (velocity dispersion 
larger than the virial velocity), which are usually not 
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included in earlier studies. They find a correlation 
consistent with a = 0.4, and a large total scatter of 
almost one order of magnitude in line width. The 
value a ~ 0.4 for the exponent of the line width size 
relation corresponds to a power spectrum of turbu- 
lence oc k~^'^. 

Miesch & Bally (1994) have estimated the power 
spectrum of turbulence in molecular clouds by com- 
puting the autocorrelation and structure functions of 
emission line centroid velocities. Their results cor- 
respond to an average exponent for the line width- 
velocity relation a = 0.43, or a power spectrum with 
(5 = 1.86. Previous attempts to measure the power 
spectrum of turbulence in molecular clouds with the 
same method had provided much shallower spectra 
(Kleiner & Dickman 1987; Hobson 1992). 

A new method to estimate the power spectrum of 
turbulence in molecular clouds has also been recently 
proposed by Brunt & Heyer (2002), using the Princi- 
pal Component Analysis by Heyer & Schloerb (1997). 
The method has already been applied to 23 molecular 
clouds in the outer Galaxy, and the result is a power 
spectrum with exponent varying from cloud to cloud, 
in the range 1.72 < j3 < 2.9, with a typical error of 
0.08. The average exponent is /3 = 2.17 ± 0.31. This 
method has been calibrated using stochastic fields 
(Stutzki et al. 1998), with purely random phases and 
no correlation between density and velocity. Correla- 
tions in real turbulent flows are likely to be important, 
and could affect the calibration of this method. 

4. The PDF of Mass Density 

The study of Probability Density Functions (PDF) 

in turbulent flows has received increasing attention 
over the last few years. PDFs can provide impor- 
tant information complementary to power spectra. A 
well known example of a combined use of power spec- 
trum and PDF, limited to linear density fluctuations, 
is the Press-Schechter model of the galaxy mass dis- 
tribution (Press & Schechter 1974). 

A number of numerical studies have established 
that the PDF of mass density in isothermal turbulent 
flows is well approximated by a Log Normal distri- 
bution (Vazquez- Semadeni 1994; Padoan, Nordlund 
& Jones 1997; Scalo et al. 1998; Passot & Vazquez- 
Semadeni 1998; Nordlund & Padoan 1999; Ostriker, 
Gammie & Stone 1999), which can be understood an- 
alytically (Nordlund & Padoan 1999). A highly ra- 
diative turbulent flow develops a complex system of 



interac;ting shocks that are able to fragment the mass 
distribution into a random network of dense cores, 
fllaments and sheets and low density "voids" , with a 
large density contrast. The intermittent nature of the 
Log-Normal PDF of mass density means that most of 
the mass concentrates in a small fraction of the total 
volume of the simulation. 

The Log-Normal distribution may be written as: 
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where n' is the number density in units of the average 
density no, 

n' = n/no, (7) 



the mean of the logarithm of density. Inn', is deter- 
mined by the standard deviation of the logarithm of 
density, a: 



Inn' = — ■ 
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which is found to be a function of the rms Mach num- 
ber of the flow A4 : 



or, for the linear density: 



CTp = bM 



(9) 
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where b « 0.5, from numerical experiments. The 
standard deviation of the linear density distribution 
thus grows linearly with the rms Mach number of the 
flow (Nordlund & Padoan 1999; Ostriker, Gammie & 
Stone 1999). 

5. The Mass Distribution of Dense Cores 

In this section we derive the mass distribution of 
dense cores based on the two following assumptions: 

i) The power spectrum of turbulence is a power law; 

ii) The typical size of a dense core scales as the thick- 
ness of the postshock gas. 

If the typical size of a dense core is comparable to 
the thickness of the postshock gas, A (§ 2), its mass 
m is: 
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where we have used the jump conditions (|]) and (^). 
Equation (|l|) shows that on a given scale L, the mass 
of dense cores is proportional to the total mass avail- 
able {pqL^) divided by the second power of the Mach 
number on that scale (Ml). 

Given the power spectrum (j^), the rms velocity, 
on the scale L is 

cr« cx L", (12) 



where 



(3-1 
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The typical shock velocity on the scale L is therefore 
<7v{L), and the shock Mach number is given by @), 
where v is replaced by ay{L). Substituting this scale 
dependent expression of AIq into ( [TTI ) one obtains 
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(14) 



where Lq is the (large) scale where the turbulent ve- 
locity is Vq and the rms Mach number is A4a,o- 

The smallest scale where significant density fluc- 
tuations may be expected is approximately the scale 
where the Mach number is of order unity: 
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In MCs with a mass Mq = poLg = 10'' Mq a typical 
value of the Mach number is Ma,o ~ 10- We then 
obtain to™„ ~ 0.0003 Mq, for (3 = 1.8. 

If Lq in equation (^4|) is defined as the largest scale 
of the turbulent flow (the scale of turbulent energy 
injection), and we take L = Lq, we obtain an estimate 
of the mass of the largest cores formed by turbulent 
fragmentation. 
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where A4a,o is the rms Mach number on the largest 
turbulent scale. In MCs with a mass Mq — PqLq = 
Mq and Mach number Ma,Q = 10, rn„iax ^ 
100 Mq. 

To arrive at the distribution of core masses we first 
consider the mass distribution of a completely self- 
similar case and then consider the modification arising 
from the Mach number dependence. 

The completely self-similar case leads to equal 
mass contributions from each logarithmic interval 



(IMF slope -1), as shown by Elmegreen (1997), who 
also demonstrated that this result does not depend on 
the fractal dimension of the self-similar distribution. 
To illustrate this result, and to derive the modifica- 
tion of it due to the Mach number dependence, it is 
useful to consider the setup and interpretation of nu- 
merical experiments. 

Numerical turbulence experiments are in a certain 
sense scale-free; the density is usually rescaled so that 
the average density (p) = 1 and the size of the box is 
scaled so that L = 1. The only distinguishing, scale 
dependent properties that remain after such scalings 
are the (sonic and Alfvenic) Mach numbers; from Lar- 
son's relations (Larson 1981) we expect larger scales 
to correspond to larger rms Mach numbers. 

To recover the completely self-similar result we 
consider two experiments with identical initial con- 
ditions (also with respect to Mach numbers and av- 
erage gas density), but interpreted at different scales 
Li, and L2 > Li. The total mass in the "large scale" 
experiment is obviously (^2/^1)^ larger than in the 
"small scale" experiment. The cores in the large scale 
experiment would be equal in number, but heavier by 
the ratio {L2/L1Y than cores in the small scale exper- 
iment. On the other hand, the total number of cores 
in the small scale experiment is (^2/^1)'^ larger than 
in the large scale experiment, if the same total mass 
is used in the two cases (that is if cores from a num- 
ber (LilLif' of small scale experiments are counted 
together). The result is therefore a total number of 
cores that depends on scale as: 



(17) 



in agreement with Elmegreen's (1997) result. 

When the Mach number dependence on scale is 
taken into account the result is that the larger scales 
contribute relatively less massive cores, because of the 
scaling relation (|lj). We assume that the number of 
cores per scale L still scales as . Combining the 
relations (|lj) and ([l^) we obtain: 



iV(m)dlnmcxm '^^^^ '^Mlnm. 



(18) 



If the spectral index is consistent with the observed 
velocity dispersion-size Larson relation (Larson 1981) 
and with our numerical and analytical results (Boldyrev, 
Nordlund & Padoan 2001), then (3 = 1.7 A and the 
mass distribution is 



N{m)A In TO oc m ^'^'^d In to. 



(19) 
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which is almost identical to the Salpeter stellar IMF 
(Salpeter 1955). 

6. The Mass Distribution of Collapsing Cores 

The mass distribution of dense cores has been com- 
puted assuming that the pre-shock density is uq and 
the postshock density MaTio, where Ma is scale de- 
pendent. A more precise computation should include 
the effect of the probability distribution of the value of 
Ma Ski each scale, or the overall effect of the statistics 
of the turbulent velocity field, which is the genera- 
tion of a Log Normal PDF of mass density (sec § 4). 
This is necessary to compute the fraction of dense 
cores that are gravitationally unstable and collapse 
into protostars, since dense cores can be significantly 
denser than their average density predicted by the 
scaling laws. While most of the large cores will be 
dense enough to collapse, the probability that small 
cores are dense enough to collapse is determined by 
the PDF of mass density. Because of the intermit- 
tent nature of the Log Normal PDF. even very small 
(sub-stellar) cores have a finite chance to be dense 
enough to collapse. 

We write the thermal Jeans' mass as: 
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is the Jeans' mass at the mean density tiq. The distri- 
bution of the Jeans' mass is obtained from the PDF of 
density assuming constant temperature as in Padoan, 
Nordlund & Jones (1997): 
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where mj is in solar masses, and: 



A = ln m J — In n' 
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The fraction of cores of mass m with gravitational 
energy in excess of their thermal energy is given by the 



integral of p{raj) from to m. The mass distribution 
of collapsing cores is therefore 



7V(m)dlnmocm-^/(^- 



/ p{mj)dmj rflnm 
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(24) 

The mass distribution is plotted in Figure 1, for 
(3 = 1.8. In the top panel the mass distribution is 
computed for throe different values of the largest tur- 
bulent scale Lo, assuming Larson type relations (Lar- 
son 1981) to rescale no and Ma,o according to the 
value of Lq. The mass distribution is a power law, 
determined by the power spectrum of turbulence, for 
masses larger than approximately 1 m©. At smaller 
masses the mass distribution flattens, reaches a maxi- 
mum at a fraction of a solar mass, and then decreases 
with decreasing stellar mass. Collapsing sub-stellar 
masses are found, thanks to the intermittent density 
distribution in the turbulent flow. The middle and 
bottom panel of Figure 1 show the dependence of the 
mass distribution on the rms Mach number of the flow 
and on the average gas density respectively. 

The magnetic critical mass is derived in the next 
section. We have not used it here to obtain the mass 
distribution of collapsing cores because the thermal 
Jeans's mass is a more strict condition for collapse. 
The magnetic critical mass depends on the core mor- 
phology in relation to the field geometry and on the 
magnetic field strength that correlates with the gas 
density with a very large scatter (see below). It is pos- 
sible therefore that magnetic pressure support against 
the gravitational collapse affects the shape of the mass 
distribution, but only as a secondary effect. 

7. The Stellar IMF 

Observations show that the stellar IMF is a power 
law above 1-2 m©, with exponent around the Salpeter 
value X = 1.35, roughly independent of environment 
(Elmegreen 1998, 2000), gradually flattens at smaller 
masses, and peaks at approximately 0.2-0.6 m© (Hil- 
lenbrand 1997; Bouvier et al. 1998; Luhman 1999; 
Luhman & Rieke 1999; Luhman 2000; Luhman et al. 
2000). The shape of the IMF below 1-2 mg, and par- 
ticularly the relative abundance of brown dwarfs, may 
depend on the physical environment (Luhman 2000). 

The scalings discussed above result in a mass distri- 
bution of dense cores consistent with the stellar IMF 
for masses larger than 1 m©, without invoking a sam- 
pling rate proportional to the free fall time, or "com- 
petition for mass" as in Elmegreen (1997, 1999). Two 
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Fig. 1. — Mass distributions of gravitationally un- 
stable cores from equation (|2^). Top panel: Mass 
distribution for different values of the largest turbu- 
lent scale Lq, assuming Larson type relations (for 
rescaling no and A^a.o with Lq), Tq = 10 K and 
/3 = 1.8. Middle panel: Mass distribution for dif- 
ferent values of A^a.Oj assuming no = 500 cm~'^, 
To = 10 K and /3 — 1.8. Bottom panel: Mass distri- 
bution for different values of no, assuming AIq O — 10; 
To = 10 K and /3 = 1.8. The mass distribution peaks 
at approximately 0.4 m©, for the values M.afl = 10, 
no = 500 cm-3, To = 10 K and /3 1.8, typical of 
nearby molecular clouds. 



conclusions are possible; either there are effects in ad- 
dition to those considered by Elmegreen, and they 
all happen to cancel each other, or else additional 
effects are not important in the first place. In the 
spirit of Occam, let's consider the latter possibility. If, 
as argued elsewhere by Elmegreen (2000a), star for- 
mation essentially happens in a crossing time, then 
we may indeed see only one generation of stars be- 
ing produced at each scale, rather than the repeated 
process implied by scaling with the local dynamical 
time. The picture thus is one where a particular MC 
forms as a consequence of the random intersection of 
counter-streaming, super-sonic motions (Ballesteros- 
Paredes, Hartmann & Vazquez-Semadeni 1999; Hart- 
mann, Ballesteros-Paredes & Bergin 2001), internal 
turbulence creates the distribution of core masses de- 
rived above, and the cores are then grabbed by grav- 
itation to form one generation of stars. Energy feed- 
back from stars subsequently disperses the cloud be- 
fore the process has time to repeat. 

In the process envisaged above, turbulent fragmen- 
tation is responsible for creating the core mass dis- 
tribution, while gravity is only responsible for the 
collapse of each protostar. The flattening and the 
turn around of the IMF is also easily accounted for 
in such a model. While scale-free turbulence gen- 
erates a power law mass distribution down to very 
small masses, only cores with a gravitational binding 
energy in excess of their magnetic and thermal en- 
ergy can collapse. The shape of the stellar IMF is 
then determined by the PDF of gas density, that is 
by the probability of small cores to be dense enough 
to collapse. The mass distribution of collapsing cores 
derived in the previous section and based on the Log- 
Normal PDF of mass density is indeed consistent with 
the observed IMF. 

The scaling of the mass where the IMF peaks can 
be derived without a knowledge of the PDF of mass 
density, using the scaling laws and the definition of 
the critical mass for collapse. We first consider the 
magnetic critical mass, 

™B=mB,o(^) (^) ' (25) 

where tob^o is the magnetic critical mass at the aver- 
age number density no. 




(McKee et al. 1993). Padoan & Nordlund (1999) 
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have shown that super-sonic and super- Alfvenic tur- 
bulence generates a correlation between gas density 
and magnetic field strength, consistent with the ob- 
servational data. The two most important properties 
of such a B-n relation are the very large scatter, and 
the power law upper envelope {B cx n^^'^). More re- 
cently, Padoan et al. (2001b) have computed the mag- 
netic field strength in dense cores produced in numer- 
ical simulations of self-gravitating, super-sonic and 
super-Alfvenic turbulence. They found typical field 
strength as a function of column density in agree- 
ment with new compilations of observational samples 
by Crutcher (1999) and Bourke et al. (2000). Here we 
adopt the following empirical B-n relation consistent 
with our previous works: 



B = B,{^ 
Po 



0.5 



(27) 



where the exponent is 0.5, and not 0.4 as reported 
above, because we now refer to the average values of 
B inside bins of n, and not to the upper envelope of 
the B-n relation, as above. The slight steepening is 
due to the fact that the lower envelope of the B-n 
relation is steeper than the upper envelope. 

Note that (|27| ) is a statistical relation between the 
average magnetic field strength as a function of core 
density, while Eqs. (|l|) and (||) are fundamental shock 
jump relations. Flows along magnetic field lines are 
able to alter the ratio of mass to magnetic flux in 
regions upstream of shocks, resulting in the ensemble 
that (^^ characterizes. Both relations are valid, and 
enter into the derivation of the critical mass. 

We find the critical mass by imposing m — niB, 
where m is given by equation (|lj), and ms by equa- 
tion (Est), 



(28) 
(29) 



^t3\ (/3-1)/(15-3/3) 

For P = 1.8, we get 

mB,c ~ msfiMao^^^- 



The critical mass is therefore typically a few times 
smaller than the critical mass at the average density. 
For Ma,o = 10, no = 10^ cm~^ and Bq = 8 ^G, 
the critical mass is approximately 2 m©.^ The prob- 

^This value of Bq = 8 fiG provides a normalization of the B—n 
relation consistent with the results of our super-Alfvenic nu- 
merical simulations discussed in Padoan & Nordlund (1999) 
and in Padoan et al. (2001b). 



ability that cores smaller than this mass are larger 
than their critical mass decreases with decreasing core 
mass, which could produce some flattening of the stel- 
lar IMF. Because of the large scatter in the B-n rela- 
tion the magnetic critical mass does not deflne a sharp 
cut-off in the IMF, but rather a gradual flattening. A 
sharper mass scale is defined by the Jeans' mass in 
the cores, 

(30) 



TOj,c ~ nijfiM^l^''^ '''' 



where toj.o is the Jeans' mass at the average density 
defined in For Ma.o = 10, = 500 cm^^. 

To = 10 K and /3 = 1.8, the thermal critical mass 
is approximately 0.4 m©. Therefore, using physical 
parameters typical of nearby molecular clouds, the 
present model predicts that the IMF should gradu- 
ally flatten below approximately 2 m© and peak in 
a (logarithmic) neighborhood of 0.4 m©, due to in- 
creasing thermal pressure support at smaller masses 
(cores smaller than their Jeans' mass are not included 
in the mass distribution). 

This result is consistent with the analytic expres- 
sion of the IMF derived in the previous section based 
of the density PDF. As can be seen in Figure 1, the 
IMF peaks at approximately 0.4 m©, for Ma.o — 10, 
no = 500 cm-3. To = 10 K and /3 = 1.8. 

7.1. The Largest Stellar Mass 

In § 5 we have estimated the largest mass of dense 
cores formed by the process of turbulent fragmenta- 
tion, given by the expression (|l6|). The largest stel- 
lar mass should also be of the order of nri^ax- As- 
suming that the turbulent velocities are of the order 
of the virial velocities in the parent molecular cloud, 
and adopting the Larson relation po c>c L^"^ (Larson 
1981), we obtain nimax oc M^i^^^ ^. This is close to 
the empirical relation nimax oc M°j *^^ (Larson 1982), 
which is found for a cloud sample known to follow the 
above size-density relation. It is likely that the true 
exponent of this relation is slightly larger than the 
value found by Larson (1982), because the lifetime of 
the most massive stars is comparable to, or shorter 
than, the lifetime of their parent molecular clouds, 
and therefore the probability of observing the most 
massive stars decreases with increasing mass. 



^According to the Larson relations larger clouds have lower av- 
erage gas density than smaller clouds. This does not mean that 
the most massive stars are formed in lower density regions. In 
each cloud stars of different mass are assumed to form, on av- 
erage, from pro-shock gas at the mean cloud density. 
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Elmegreen (1993, 1997) has argued that the largest 
stellar mass is related to the mass of the parent cloud 
for purely statistical reasons: the larger the cloud 
mass is, the higher the probability of populating the 
high mass tail of the IMF. Such a statistical argument 
is correct only if the normalization of the IMF is inde- 
pendent of the total cloud mass. The scalings derived 
above show that this is not the case; larger clouds 
in general have larger velocities and therefore form 
relatively fewer stars of a given mass. This is also 
indicated by observations, since it is commonly found 
that larger clouds have lower star formation cfhcicncy 
than smaller clouds. However, since the star forma- 
tion efficiency also depends on the low mass cut-off of 
the IMF, and since age differences also may enter, it 
would have been hard to draw firm conclusions from 
observations alone. 

8. Discussion 

In Section ^ we found that the mass distribution 
of dense cores formed by turbulent fragmentation has 
a power law shape, with a slope consistent with the 
stellar IMF at intermediate and large masses, and in 
Section ^ we found the distribution of the mass of 
collapsing cores, which has a shape consistent with 
the observed IMF also at small masses. Three as- 
sumptions were essential for the derivation, and we 
would here like to emphasize what these assumptions 
are, and point out that they are independent of each 
other. 

The first assumption is that approximate selfsimi- 
larity holds for supersonic turbulence of varying Mach 
number, to the extent that the scaling of the number 
of cores with upstream scale L scales as , also 
when the Mach number varies with scale L. This 
is a reasonable assumption because the compressed 
cores occupy only a tiny fraction of the volume where 
their mass was collected from. The second assump- 
tion is that the mass of the compressed cores scales 
as the cube of the thickness of the shocked gas, which 
again scales as the inverse of the Mach number. This 
assumption is independent of the particular density 
and scale of the upstream flow, and rests only on the 
fundamental shock jump condition and on the core 
shape being independent of Mach number. The third 
assumption is that the Mach number of the upstream 
gas depends on the scale L of the upstream flow as 
L", where a is essentially the index in Larson's (1979, 
1981) velocity-size relation. 



It is remarkable that, under these assumptions, an 
IMF slope consistent with Salpeter's (1955) value is 
obtained as a direct result of the observed velocity 
dispersion-size relation of ISM turbulence and the 
fundamental jump conditions for isothermal MHD 
shocks. Provided that some general criteria are met, 
the slope of the stellar IMF is thus independent of 
the physical conditions in the star-forming clouds, as 
also indicated by the observations (Elmegreen 1998, 
2000). 

We have also interpreted the gradual flattening of 
the IMF around 1-2 Mq and its peak at approxi- 
mately 0.3-0.5 Mq as the effect of thermal (possibly 
also magnetic) support against the gravitational col- 
lapse. The mass distribution of collapsing cores, com- 
puted on the basis of the PDF of mass density and the 
thermal Jeans' mass, is consistent with the observed 
stellar IMF, down to sub-stellar masses. 

We have shown that the IMF at low stellar masses 
is sensitive to the average physical properties of the 
star forming gas, such as the rms Mach number and 
the average gas density, also consistent with the ob- 
servations (Luhman 2000). An even larger number 
of sub-stellar objects would be predicted if a signif- 
icant fraction of them originate as close companions 
of more massive stars and later on separate as in- 
dividual brown dwarfs. In the present work we do 
not study the evolution of density fluctuations smaller 
than their Jeans' mass within the collapsing back- 
ground of an unstable larger core, and therefore this 
mechanism for the formation of giant planets and 
brown dwarfs is not discussed here. 

The assumption behind the result is that the up- 
stream density is sampled from a Log-Normal density 
distribution, and that the distribution of the den- 
sities of cores therefore also is approximately Log- 
Normal. An exact Log-Normal density distribution 
requires exactly isothermal conditions, but even for 
non-isothermal conditions the central part of the den- 
sity PDF is still approximately Log-Normal (Scalo et. 
al 1998, Nordlund & Padoan 1999). 

Given the fact that turbulent fragmentation is un- 
avoidable in super-sonic turbulence, and given the 
success of the present model in predicting the correct 
slope of the stellar IMF, it is difficult to argue that 
super-sonic turbulence does not play a dominant role 
in the generation of the stellar IMF. Other processes 
such as gravitational fragmentation (Larson 1973; 
Elmegreen & Mathieu 1983; Zinnecker 1984), opac- 
ity limited fragmentation (Hoyle 1953; Gaustad 1963; 
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Yoneyama 1972; Suchkov & Shchekinov 1976; Low 
& Lynden-Bell 1976; Rees 1976; Yoshii & Saio 1985; 
Silk 1977a, b), protostar interactions and coagulation 
(Nakano 1966; Amy & Weissman 1973; Silk & Taka- 
hashi 1979; Bastien 1981; Yoshii & Saio 1985; Lejeune 
& Bastien 1986; Allen & Bastien 1995, 1996; Price 
& Podsiadlowski 1995; Murray & Lin 1996), stellar 
winds and outflows (Silk 1995; Nakano, Hasegawa & 
Norman 1995; Adams & Fatuzzo 1996), competitive 
accretion (Larson 1978; Tohline 1980; Bonnell et al. 
1997; Myers 2000) must be relatively unimportant. 

This conclusion is supported by a recent computa- 
tion of the mass distribution of dense self-gravitating 
cores in numerical simulations of self-gravitating super- 
sonic MHD turbulence (Padoan et al. 2001b). The 
result is a power law mass distribution consistent 
with the stellar IMF. Mass distributions of cores 
from numerical simulations of super-sonic turbulence, 
roughly consistent with the stellar IMF, are also re- 
ported by Klessen (2000). 

It should be noted that even the formation of 
molecular clouds, a process that some authors as- 
cribe to the random collisions of gas streams (e.g. 
Ballesteros-Paredes, Hartmann & Vazquez-Semadeni 
1999; Hartmann, Ballesteros-Paredes & Bergin 2001), 
may be a manifestation of a process analogous to the 
one discussed above, just operating on a larger scale. 
If motions on larger scales are also characterized by 
power laws, one would then expect the distribution of 
MC masses to obey a Salpeter-like scaling. However, 
such scaling would be rather difficult to probe obser- 
vationally, due to the fact that MCs are dispersed by 
the process of star formation before they are able to 
collapse as a whole, which prevents them from being 
well defined as individual objects. 

For almost twenty years estimates of the mass dis- 
tribution of dense cores in molecular clouds, based 
on molecular-line studies, found a shallow power law 
mass distribution with a single exponent in the range 
< a; < 0.70, with a typical value x = 0.5. The 
main reasons for the shallow mass distribution ob- 
tained in these works are i) the relatively low den- 
sity traced by the molecular emission lines normally 



5 (Myers, Linke & Benson 1983; Casoli, Combes & Gerin 1984; 
Blitz 1987; Carr 1987; Loren 1989; Stutzkie & Glisten 1990; 
Lada, Bally & Stark 1991; Nozawa et al. 1991; Tatematsu et 
al. 1993; Langer, Wilson & Anderson 1993; Williams & Blitz 
1993; Blitz 1993; Williams, De Geus & Blitz 1994; Williams, 
Blitz & Stark 1995; Dobashi, Bernard & Fukui 1996; Onishi et 
al. 1996; Yonekura et al. 1997; Kawamura et al. 1998). 



used (^ 10^ cm~^ for ^^CO, and ~ 10'' cm~^ for 
C^^O); ii) the limited density range probed by each 
molecular emission line; iii) the use of inappropriate 
"clump-find" algorithms; iv) the relatively low res- 
olution that allows only the selection of large cores, 
with typical mass ranging from sa 10 m© to hundreds 
or thousands uiq. Clearly, cores selected in this way 
cannot be identified with single protostellar cores, as 
discussed in Padoan (1995), where it was predicted 
that the exponent of the intrinsic mass distribution 
of protostellar cores should have been x > 1, based 
on the stellar IMF. Recently, Onishi et al. (1999) 
have obtained a sample of dense cores in the Tau- 
rus molecular cloud complex using a higher density 
tracer, H^'^CO"'', which probes a density of approxi- 
mately n = 10^ cm~'^. They found a power law mass 
distribution with exponent x = 1.5 ±0.3, in the mass 
range between 3.5 mQ and 25 m©. While most pre- 
vious determinations of core mass distributions using 
molecular-line maps are affected by the arbitrary def- 
inition of an individual core, which is far from trivial 
in the hierarchical cloud structure, the dense H^'^CO"'" 
cores found by Onishi et al. (1999) in Taurus are all 
isolated and therefore unambiguously defined. 

Recent dust continuum emission surveys, which 
also probe relatively high densities (n = 10^-10^ cm~^), 
have provided more support to the idea that the stel- 
lar IMF reflects the mass distribution of dense cores. 
The mass distribution of dense cores in p Ophiuchi 
(Motte et al. 1998; Johnstone et al. 2000), in the Ser- 
pens core (Testi & Sargent 1998) and in Orion (Motte 
et al. 2001; Johnstone et al. 2001) are found to be 
consistent with the stellar IMF. Motte et al. (1998), 
for example, obtained a power law mass distribution 
with exponent x = 1.5 (in logarithmic units such 
that Salpeter's exponent is a; = 1.35) in the range 
of masses between 0.5 m© and 3 m©, and x = 0.5 in 
the range between 0.1 m© and 0.5 m©. The results by 
Motte et al. (1998, 2001) are confirmed by Johnstone 
et al. (2000, 2001). Testi & Sargent (1998) found 
X = 1.1, in the range of masses between 0.5 m© and 
30 m0. 

Detailed studies of pre-stellar cores selected in this 
way could provide accurate volume density measure- 
ments and therefore a test for the density-size relation 
implied by our equation ([Tl|). Note, however, that 
because all quantities involved have wide distribution 
functions it is important to distinguish between av- 
erage or typical values, extreme values, and observed 
values. The latter may often be more representative 
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of extreme values than of average values. One ex- 
ample is Larson's density-size relation, which by its 
nature tends to refer to the densest cores of a given 
size rather than the average cores of that size; many 
low density cores may exist that are unobservable, or 
are ignored. 

In the context of our derivation and assumptions 
one may understand Larson's density-size relation as 
the result of combining Eq. ^ and ^ and then taking 
the extreme with respect to L; i.e., the relation rep- 
resents those extreme cases where roughly the same 
initial density gas has been compressed by shocks, 
perhaps repeatedly, to various degrees. If each com- 
pression is essentially one-dimensional one retains the 
approximately inverse relation between core size and 
core density. 

9. Conclusions 

In conclusion, we have related the stellar IMF to 
the mass distribution of dense cores formed by the 
process of turbulent fragmentation, assuming that 
only cores with gravitational energy in excess of their 
magnetic and thermal energy can collapse as proto- 
stars. Most sub-critical cores disperse back into the 
turbulent flow, and are therefore irrelevant for the 
process of star formation, as we have argued in other 
recent works (Padoan et al. 2001a, b). Previous the- 
ories of star formation (sec Shu, Adams & Lizano 
1987) assume that stars of small and intermediate 
mass are formed from sub-critical cores. In such the- 
ories, sub-critical cores are in static equilibrium and 
evolve quasi-statically, on the time-scale of ambipolar 
drift. We have argued that such a scenario is incon- 
sistent with the turbulent nature of MCs (Padoan et 
al. 2001a, b). 

We have derived the mass distribution of dense 
cores generated by the process of turbulent fragmen- 
tation, and have found a power law mass distribution 
consistent with the Salpeter stellar IMF, for stellar 
masses larger than 1-2 uiq . This result is a direct con- 
sequence of fundamental physical properties of super- 
sonic turbulence in MCs, such as the power spectrum 
of turbulence and the jump conditions for isothermal 
MHD shocks. We have also shown that another fun- 
damental physical property of turbulent flows, namely 
the PDF of the gas density, can explain the shape of 
the IMF at smaller stellar masses, and the formation 
of gravitationally unstable cores of sub-stellar mass. 

The main results of this work are: i) The power 



law stellar IMF at masses larger than 1-2 mQ is the 
result of the near self-similar nature of inertial range 
super-sonic turbulence; ii) The low mass roll-over and 
cut-off of the IMF is caused by the combined effects 
of the thermal support of the smallest cores against 
gravitational collapse and the density PDF of super- 
sonic turbulence; iii) The mass-scaling of the peak 
of the IMF may be expressed as a function of the 
physical parameters of turbulent star-forming clouds; 
primarily their rms turbulent velocity, temperature, 
magnetic field strength and average density; iv) For 
physical parameters typical of nearby MCs the peak 
of the IMF is predicted to be around 0.3-0.5 m0; v) 
The slope of the IMF at masses larger than 1-2 uiq 
is determined by the inertial range spectral index of 
super-sonic turbulence and the jump conditions for 
isothermal MHD shocks; vi) For magnetically domi- 
nated jump conditions, which are applicable to typical 
molecular cloud conditions, the IMF spectral index is 
equal to 3/(4 — /3), where /3 is the inertial range spec- 
tral index of super-sonic turbulence; vii) For a value 
of the inertial range spectral index /? — 1.74, consis- 
tent with Larson's relations and with new numerical 
results, the IMF spectral index is a; = 1.33, almost 
identical to Salpeter's slope. 

We conclude that the process of turbulent frag- 
mentation is essential to the origin of the stellar IMF, 
in support of the thesis that star formation can be 
viewed as the main consequence of the dissipation of 
super-sonic turbulence in molecular clouds. 
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